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We present a quantum analysis of the massless excitations in graphene with a charge impurity. 
When the effective charge exceeds a certain critical value, the spectrum is quantized and is un- 
i-^ ' bounded from below. The corresponding eigenstates are square-integrable at infinity and have a 

ly-j ' rapidly oscillatory behaviour in the short distance, which can be interpreted as a fall to the centre. 

(\J , Using a cutoff regularization, we show that the effective Coulomb interaction strength is driven to 

its critical value under the renormalization group flow. In the subcritical region, we find bound 
states with imaginary values of the energy for certain range of the system parameters. The physical 
significance of these bound states with imaginary eigenvalues is discussed. 
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In the tight-binding approximation, the excitations of graphene near the Fermi surface are described by a two 
dimensional (2D) massless Dirac equation [JQ- The charge impurities in graphene follow the laws of 2D quantum 
electrodynamics 3, 4, 5, 6] and vortex fields lead to fermion fractionalization 0, Recent fabrication of graphene 
' monolayers Q provides an opportunity to experimentally study the effects of charge impurities [lol 11 , . The 



massless Dirac excitations normally do not form bound states with the impurity. Numerical analysis of the coulomb 
impurity problem in graphene indicates the existence of bound states in the supercritical region [loj . where the 
■rj" ■ Coulomb potential strength exceeds a critical value. Semiclassical analysis also suggests that an infinite family of 
[ quasibound states appear in the supercritical region The Fermi velocity vp in graphene is approximately 10^ 

• 2 

lO ' m/s and the effective fine structure constant ag = ~ 2.5. This corresponds to the critical value of the effective 
, impurity charge ~ 1. Measurement of local density of states using scanning tunneling spectroscopy (STS) has recently 
led to experimental evidence for quasibound states in this system 

In this Letter we shall present a full quantum mechanical analysis of the bound states of massless excitations in 
graphene in the presence of a Coulomb potential. In the supercritical regime, we shall show that there exists an 
infinite number of bound states. We establish that the spectrum is unbounded from below, which is a characteristic 
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\^ feature of the "fall to the centre" [IJ, llJ]. This feature was absent in the semiclassical analysis of [121. We also find 
5^ , that the spectrum in the supercritical region is labelled by a parameter which is not determined by the theory. The 
significance of this parameter can be understood as follows. The two dimensional Dirac description of graphene is a 
long wavelength approximation of the underlying dynamics. In the supercritical region, the power law behaviour is 
likely to fail as the short distance effects become important @ . Our result suggests that the effect of the breakdown 
of the power law potential at short distance can be parametrized by a single quantity. This parameter is not predicted 
by the theory itself, but must be determined empirically. The ratio of the eigenvalues is however independent of this 
parameter. 

A special feature of graphene with charge impurities is that the nonlinear screening effects drives the effective charge 
to the critical value [llj . In our formalism this is also present and can be interpreted as a renormalization group flow, 
which is well known for singular inverse square interactions ■ We will show that the corresponding beta function 
has an ultraviolet stable fixed point at the critical value of the Coulomb potential. 

For the massless two dimensional Dirac equation, there is no bound state in the subcritical region. This statement 
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is true only for real value of the energy. We find that for certain range of the parameters, the system admits square- 
integrable bound states with eigenvalues ±i. Using von Neumann's theory of self-adjoint extensions [l6j . we show 
that these states lead to a one parameter family of boundary conditions for the relevant range of system parameters. 
These generalized boundary conditions are expected to affect the phase shifts and the S-matrix in the scattering 
sector. This situation is analogous to the occurrence of anomalies for singular potentials 17, 1^ . 

The 2D massless Dirac equation in the presence of a Coulomb potential V{p) = can be written as 



hvf 





-idx + dy 



-idx - dy 




X=[E- V{p)]x, 



(1) 



where p is the radial coordinate on the two dimensional x-y plane. Let (f) denote the corresponding polar angle on the 
plane. Following JTj, the separation of variables can be performed using the ansatz 



xip, 0) 
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where 



m is the half integer azimuthal quantum number and 



E 



Khvp ' 



(3) 



with n as the effective dielectric constant. The minimum magnitude of the azimuthal quantum number is m = i and 
from ^ we see that the corresponding critical value of /3 is given hy j3c — \- Assuming that graphene electrons are 
the only source of screening gives k « 5 fs!], which together with Pc = \, gives the critical value of Z as Zc ^ 1 11 1. 
This feature makes the system interesting from an experimental point of view. 
The Dirac equation ^ with the ansatz ([2]) leads to two coupled equations. 



du r, n 1 \ 

p— — h (s + zp + 2ikp)u — mv 
dp 

dv , 
p— — h (s — ip)v — mu 
dp 



0, 
0. 



Defining z = —2ikp, from (4) and (5) we get 

d^v 



dv 



dz^ 



+ (1 + 2.S - z)— - (s - if3)v = 0, 
dz 



(4) 
(5) 

(6) 



which is of the form of Kummer's equation for the confluent hypergeometric function [10|. We now define a new 
function rj as 



e-'^P 77. 



Using ([71) in ([6]), we get 



where a = + i(3. Eqn. 



d^rj 
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^ = 0, 



(7) 



(8) 



has the form of Whittaker's equation [H 



We shall first analyze the supercritical case where /? > /3c- In this case, s = Hp, p — \/ P"^ ~ rn? and Eqn. ([8]) can 
be written as 



dz2 



z 



The general solution of ([9]) is given by 
where Ai and A2 are constants and 



M, 



a, ±1/1 



-%z^±^t^M{±ip -ip,l± 2ip, z), 



(9) 
(10) 
(11) 
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with M on the right hand side of pTjl denoting the confluent hypergeometric function p^ . 

We shall now construct solutions which are square integrable at infinity. Let us first note that for the functions 
F{p) and G(p) in the measure of integration is given by pdp. Now consider the contribution of the function 77 in 
(fTU| to the quantity J \F\'^pdp at asymptotic infinity. Using the asymptotic formula (13.5.1 of [31) for the confluent 
hypergeometric function, we find that as |z| = | — 2ikp\ 00, 

M„,±.,(z) e'^^e-(^^-^)(-2zM^+' V,7^^y^, + e~^^^ {-mp)-'^-^^ ^j^^^^. (12) 

i{l±ip + ip) l[±ip—ip) 
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Using ([21), ^ and (fT2)) . it is easy to see that as p 00, the contribution of the first term on the rhs of (fT2l) to \F 
behaves as ~ ^ while that from the second term behaves as ~ Hence, the first term in the rhs of (|12p provides 
a divergent contribution to / \F\^pdp while that of the second term is convergent. We now consider the part of the 
solution 77 in (fTO|) . arising from the first term in the rhs of (fT2|) . which leads to the divergence of the norm of F. This 
is denoted by 77div and is given by 



?7di- 



■i'^P{-2ikp)^+'f^e-~''^ 



^^^^ r(l + 2z^) Y{\-2^p) 

' r(i + ip + i/3) ^ ' r(i -ip + i/3) 



(13) 



We now choose A2 — ^le^^^ and denote r(i+t;^+1ff) ~ ^e'^'^ ■ Therefore, as p — > 00, we get 

?7div 2Aie'''P{^2ikp)^+''^e''^''-'^kcos{e + -f). (14) 
We can make Tydiv = if we impose the condition that 

cos(6' + 7) = or 6* = -7 + ^tt. + tt, n e Z. (15) 

This gives a quantization of the parameter 9. The quantization of 9 ensures that as p — > 00, J \F\^pdp will remain 
finite. Knowing 77, the components u and v of the solution to the Dirac equation can be obtained from ([7]) and (5). 
They are square-integrable and represent the bound state solution. 

We now proceed to obtain the quantized energy levels. For that purpose, consider the short distance limit of (jlOp . 
Denoting r|i[l^^;j;.^] = Ce'^'', we get, as z ^ 0, 

77 - [(1 + C) cos(6l -5 + p\nz) + {l~C) sm{9 -5 + p\nz)]. (16) 

The short distance behaviour of the system can also be inferred by looking at the corresponding indical equation. 
In this case, the roots of the indical equation are given by ^ ± ip. Using these two roots, the wavefunction at short 
distance can be obtained as 



77 ~ \/-2ip 



- ^J-2ip [(1 + C) cos(B + p ln(-2ip)) + i (1 - C) sin(B + p ln(-2ip))] , (17) 

where _B is a real constant. The short distance behaviour of the wavefunction obtained from p6p and (|17p must 
match. This is possible provided 

9 + 8^p\\\k = B. (18) 
Thus wc obtain the quantized energy eigenvalues as 

E^^-^^-"—— , 77 eZ, (19) 

where A — ^"^"^-^ is a constant. The quantity 8 depends on the system parameters p and /3. B on the other hand is 
a real constant which is not fixed by the theoretical analysis. This leads to a one parameter family of inequivalent 
spectra in the supercritical regime. As mentioned before, B encodes the effect of short distance physics, which are 
expected to be important as we approach length scales of the order of the lattice spacings. The analysis suggests 
that it is not important to know the the details of the short distance interactions and their effect on the spectrum 
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appears through a single parameter B, which should be determined empirically. The ratio of the various energy levels 
are independent of the parameter B. It may be noted that the above quantization has been carried out with only 
the function v(p). The function u{p) can be determined from (5) and it can be shown that it does not change the 
quantization condition. 

We have thus found an infinite number of bound states in graphene containing a impurity, when the effective charge 
exceeds the critical value. These states are square-integrable at infinity and have a rapidly oscillatory behaviour in the 
short distance limit, which is a feature of the fall to the centre. Our analysis is fully quantum and nonperturbative, 



valid for any value of P > Pc- In that sense we go beyond the treatment in [l^]. The semi-classical analysis in [l^] 
yields quantized energy levels only for n > whereas we obtain quantized energy levels for all integer values of n. 
Due to this, the spectrum we get is unbounded from below, which is an important feature of the supercritical region. 
We also find an accumulation point at zero energy, which is seen in the exact numerical analysis of the system as well 

Another interesting feature of graphene is that the screening in the supercritical region drives the effective charge 
to the critical value In our formulation this is achieved when /i — s- from the supercritical region. Within the 
scope of the quantum mechanics discussed here, this effect can be viewed as a renormalization group flow arising from 
the short distance effects. In the short distance limit, the eigenvalue equation ([9]) assumes the conformal form given 
by 



V = At?, (20) 



with the eigenvalue A = 0. Eqn. for a general A has been analyzed in 15| and provides an example of renor- 
malization in quantum mechanics. In the supercritical region, the spectrum is unbounded from below. In order to 
regulate this divergence, we introduce a short distance cutoff at p = a and impose the boundary condition that the 
wave-function vanishes below the cutoff. With such a boundary condition, the eigenvalues of (pO)) can be calculated 
analytically for small values of p,, leading to a finite bound state spectrum for A given by [15I | 



A. 



(21) 



where a is the Euler's constant and n = 1, 2, oo. Note that for any finite value of the cutoff, as n — + oo, A„ — > 
and the zero eigenvalue depicts an accumulation point for this system. In this limit we recover the short distance 
form of (O. The eigenvalues A„ also explicitly depend on the cutoff and they diverge as the cutoff is removed. 
This indicates that the breakdown of the 2D Dirac description of graphene at short distances @. In the spirit of 
renormalization group analysis, we now make the coupling p a function of the cutoff, i.e. take p — p{a). In order 
to find the dependence of p on the cutoff, we demand that the zero eigenvalue remains unchanged as the cutoff is 
removed. This leads to a beta function for the coupling p given by [l5| 

M = -a^ « -A^'. (22) 

This beta function has an ultraviolet stable fixed point at p ~ 0. In other words, the effective charge is driven to its 
critical value as the cutoff is removed. We therefore recover the same result obtained from the analysis of nonlinear 



effects on screening [ll|. 



We shall now analyze graphene with charge impurity in the subcritical region with (3 < [3c- s is now a positive 
real number, s — Q denoting the critical value. In what follows we shall restrict our attention to s > 0. If the 
corresponding Dirac equation ([1]) had a mass term, there would be bound states in the subcritical region, with the 
bound state energy proportional to the mass Q . In our case we do not have the mass term in ^ and would normally 
not expect any bound state, at least with real energy. For the moment, let us consider the system where the quantity 
k is not real. In particular, we ask the question whether there exist bound states with values of fc = ±j. We shall 
discuss the significance of such solutions later. 

Consider first the case with k = +i. In this case, z = —2ikp — 2p and a possible solution of ([6]) is given by 



From (5), we get that 



v+ = U{s + 2s, z) = U{s -il3,l + 2s, 2p). (23) 

(s - iB) 

^ ^ Pis + 2s, 2p) - 2pU{s -1/3+1,2 + 2s, 2p)] , (24) 



5 



where we have used the relation U'{a, b, z) = —aU{a + l,b + 1, z) [1^. Using (|23p and in we obtain 



{m + s - 1 + 2., 2p) - - z/5 + 1, 2 + 2s, 2p) 

m TO 



(25) 



Since p is positive, as p — > oo, U{a,b,2p) ^ p Therefore we see that as p ^ oo, 0. Next consider the 

short distance Hniit of (p) . Using the relation 



U{a,b,y) 



sin(7r6) 



M{a,b,y) ^_,M{l + a-b,2-b,y) 

V 



T{l + a-b)T{b) 



r(a)r(2-6) 



(26) 



(27) 



Thus, for the range < s < ^, we see that F+{p) is a square integrable function. Similar analysis shows that for 
< s < ^, the entire radial wave-function is square integrable. Thus, for fc = i, we have a single square integrable 
bound state when < s < ^. We denote the corresponding radial bound state with 

Let us now consider the case when k = — z, where z — ~2ikp — ~2p. In this case, a possible solution of ([6]) is given 



and noting that as y ^ 0, Af(a, 6, y) ^ 1 [1^, we get that as p ^ 0, 

\F^\^pdp^ / p^^^dp+ other convergent terms. 



by [IS 



From (5), we get that 



e^f7(l + s + i(3,l + 2s, -z) = e'^PUil + s + i/?, 1 + 2s, 2p). 



'^^ U{1 + s + z/3, 1 + 2s, 2p) - 2p ^^^'' '^^ U{2 + s + tP,2 + 2s, 2p) 



TO 



TO 



From these we see that 
F-ip)^e-P 



(J!l±l^pis-^.)U{1 + s + z/3, 1 + 2s, 2p) - 2(2 + ^-»/?) ^(.+i)^(^ - z/3 + 2, 2 + 2s, 2p) 

TO TO 



(28) 



(29) 



(30) 



From ([50)1 we notice that as p ^ oo, F-{p) 0. In the short distance limit, we can again show that F_ and the 
corresponding entire radial wave- function is square integrable when < s < ^. Thus, for k = —i also, we have a 
single square integrable bound state when < s < ^ . We denote the corresponding radial bound state with tp_ . 

We now comment on the physical significance of the square- integrable bound states ip± with k = ±i, which exist 
when < s < Their significance can be understood following von Neumann's treatment of self-adjoint operators 
in quantum mechanics [l6|. The original radial Dirac operator for graphene, denoted by Hp is defined with the 
boundary condition that the wave-function vanishes at the location of the charge impurity given by p = 0. This 
boundary condition defines a domain of the Dirac operator denoted by D{Hp), whose elements vanish at p — 0. The 
corresponding adjoint operator Hj^ in this case has the same differential expression as Hp. Let n± denote the number 
of linearly independent square integrable solutions of the eigenvalue equation for Hj,, with eigenvalues zti. In our 
case, we have found that 7i_|- = n_ = 1 when < s < ^ and n± — otherwise. The theory due to von Neumann 
[l6t implies that for < s < ^, the radial Dirac operator Hp is not self-adjoint in the domain D{Hp). It furthermore 
says that Hp would be self-adjoint in a different domain, denoted by D^^Hp), whose elements can be written as 
^ -\- Q{ip+ -\- e^^ip-) where ^ € D{Hp), Q is a constant and u € R (mod 27r). In effect, this provides a one parameter 
family of self-adjoint extensions, or equivalently boundary conditions, which are labelled by a single real parameter 
V g [0,27r]. Thus, the physical significance of the bound states with fc = ±i is that they lead to a more general 
class of boundary conditions which are allowed by the rules of quantum mechanics. The quantity v which classifies 
the different boundary conditions is called the self-adjoint extension parameter. It is expected that the phase shifts 
and the S-matrix in the scattering sector of the Dirac operator |T]) would depend explicitly on the parameter v for 
< s < i [H. Thus, even though the solutions with k = ±i are not physical bound states, they do play a role in 
determining the spectrum of the problem when < s < ^. 

In summary, we have presented a quantum analysis of bound states in graphene in the presence of a charge impurity. 
In the supercritical region, we find a quantized spectrum with the energy unbounded from below. The eigenstates are 
square-integrable at infinity and have rapidly oscillatory behaviour at short distances, indicating a fall to the centre. 
They are parametrized by a single real constant, which captures the effect of short distance physics. Experimental 
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evidence for the bound states in the supercritical regime was found in [12]. It is expected that the states described 
here correspond to what was observed experimentaUy. We also show the Coulomb potential is driven to its critical 
value under a renormalization group flow. Finally, for certain range of the system parameters in the subcritical region, 
we find bound states with eigenvalues given by ±i. These states determine the allowed boundary conditions for this 
system and the corresponding self-adjoint extension parameter is expected to affect the phase shifts and the S-matrix. 
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